Effects of isospin mixing in the A = 32 quintet 
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For the A = 32, T — 2 quintet we provide a unified theoretical description for three related aspects 
of isospin mixing: the necessity of more than three terms in the isobaric mass multiplet equation, 
isospin-forbidden proton decay in '^^Cl, and a correction to the allowed Fermi Z?"*" decay of ^^Ar. 
We demonstrate for the first time that all three effects observed in experiment can be traced to a 
common origin related to isospin mixing of the T = 2 states with T — 1 states. 
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Recent measurements of the lowest T = 2 states in ^^Cl 
and '^^S have made the A — 32 multiplet the most pre- 
cisely measured T ~ 2 quintet [T], [5]. In first-order per- 
turbation theory the masses in an isobaric multiplet are 
given by the isobaric mass multiplet equation (IMME): 

M{T,)^a + bT,+c{T,)^, (1) 

where = {N - Z)/2 ^. Multiplets with T > 1 may 
require terms of higher order in that enter in second- 
order perturbation theory along with isospin mixing. The 
A — 32 multiplet requires a small but non-zero higher- 
order term, proportional to either [T] or 

Isospin-forbidden proton decay, another signature of 
isospin mixing [3], has been observed from the T — 2 
state in ^^Cl to the low- lying T = 1/2 states in ^^S with 
a decay width of 20(5) eV 

Finally, the superallowed Fermi /3+ decay of "^^Ar has 
been measured recently. Superallowed 0"*" to 0"*" Fermi 
decay is measured very precisely for many nuclei and 
provides the critical data for extracting the weak mix- 
ing angle Vud for the KCM matrix [B]. The extraction of 
Vud from the data requires a small but important correc- 
tion 5c due to isospin mixing. The correction is smallest 
for nuclei near stability, typically 0.5% or less, but can 
be larger for nuclei far from stability. A correction of 2% 
was recently measured [5]. 

These three isospin-mixing effects, usually studied and 
theoretically treated independently, are derived from the 
same origin, as evidenced by the A = 32 quintet. Calcu- 
lations for energy levels, spectroscopic factors, gamma 
decay, isospin-mixing matrix elements, and one-body 
transition densities for the multiplet were obtained in the 
sd-shell model space with the USD, USDA, and USDB 
interactions 0, [8], [9]. The shell model code OXBASH 
[TU] was utilized for a full diagonalization of the Hamilto- 
nians. The USDB interaction will primarily be used for 
explanation and illustration in this paper, with reference 
to the other interactions for meaningful comparisons. 

All calculations were carried out in proton-neutron for- 
malism. The isospin-mixing interaction is taken from the 
work of Ormand and Brown [11 , where in addition to the 
Coulomb potential, charge-independence-breaking (CIB) 
and charge-symmetry-breaking (CSB) interactions were 



TABLE I: Mass excesses of T = 2 states in the A = 32 quintet. 
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added to the USD Hamiltonian. The CIB strength was 
obtained from a one-parameter fit to the experimental 
c coefficient of the T = 1 IMME and is consistent with 
the np vs pp scattering data |3]. The CSB strength was 
obtained from a one-parameter fit to the experimental b 
coefficients. 

All 0+ states for the A — 32 quintet were calculated. 
The dominant isospin of each state was determined by 
calculating overlaps with the isospin-conserving part of 
the interactions. The lowest T = 2 state was the ground 
state for ^^Si and ^^Ar, the third 0+ state for ^^P and 
•^^Cl, and the tenth (eleventh, tenth) 0+ state for '^^S using 
the USDB (USDA,USD) interaction. 

The experimental masses are given in Table I. The 
masses for [5], [TT], ^"^Cl [T], [H], and ^'^Ar [H] are 
identical to those given in Table I of [1]. We have com- 
bined the two values for ^^V, -19232.46(15) [I3]> E] 
and -19232.78(20) [TS], [TB], into a reduced value based 
on a fit to a constant. We use the recent direct mea- 
surement at the NSCL for the '^^Si mass [12]. The circles 
with error bars in the bottom panel of Fig. 1 show the dif- 
ferences in keV between the experimental masses of the 
T = 2 states and those obtained from a fit to Eq. (1). The 
crosses display the differences obtained for the T = 2 en- 
ergies as calculated with USDB. The differences for both 
experiment and theory obtained when d[e) terms pro- 
portional to T^{T^) are added are shown in the middle 
(top) panel. For the experimental data, the value of 
each fit is given in the figure, as well as the best-fit d and 
e parameters. 

In Fig. 1, it can be seen that a much better fit oc- 
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FIG. 1: Accuracy of the IMME in the A = 32 quintet with 
three terms (bottom panel) and with an additional cubic 
(middle) or quartic (top) term. Circles correspond to experi- 
mental data, with error bars from Table I. Crosses correspond 
to theoretical calculations with the USDB interaction. Note 
the reduction in scale for the middle panel. 



curs, for both USDB and experiment, when a d coefh- 
cient is used. We repeated the procedure with Mf.xp = 
-24080.92(5) for ^^Si, combining the two indirect ^^Si 
masses [2U], [T] from Table I of p[| into a reduced value. 
While the parameters of the fit and the mass differences 
change, the conclusions are identical. Again, a d coeffi- 
cient is necessary for a reasonable fit to data, producing 
xLp = 0.58 with = 0.53(11). 

The most significant difference between theory and ex- 
periment in the bottom panel of Fig. 1 corresponds to 
the quality of the fit for ^^Ar {T^ — —2). A reduction in 
the error bar of ^^^Ar, at least to the level of '^^Cl, would 
better constrain the fit and therefore the d parameter. If 
we exclude the = —2 point, the d term can be solved 
algebraically to give dexp = 0.95(37), or d^xp — 0.41(33) 
with the indirect '^^Si mass. The evidence strongly sug- 
gests that the three-term IMME does not fit the data or 
the USDB calculations. The uncertainty in the theoreti- 
cal calculations can be assessed by comparing the results 
for the three different interactions. The USD interaction 
result patterns the USDB behavior but with larger de- 
viations than seen in Fig. 1, resulting in a greater value 
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FIG. 2: 0+ levels in ^^S given by T vs. E^. For ^^P and 
^"^Cl, the T = levels can be ignored to give an approximate 
distribution of the 0"'' states. 



of the necessary coefficient dth = 0.39. The USD A in- 
teraction cannot be corrected solely by a d coefficient, 
as large shifts in both '^^Cl and ^^S occur due to isospin 
mixing. Both the sign and the magnitude of the neces- 
sary coefficient give information about the shifts of the 
T — 2 states, which can be determined theoretically. 

Fig. 2 shows the first twenty 0+ levels in ■^^S with the 
USDB interaction, categorized by their values of T. The 
sum of the ground state energy and the excitation energy 
of the T = 2 state gives the energy to be used in the fit for 
•^^S. As seen from the graph, the T — 2 state has nearby 
states that repel it, shifting its energy. The shift of the 
T = 2 state, labeled by "a" , is approximately given by 
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where V denotes the interaction that causes isospin mix- 
ing. Ex, the excitation energy of the state, is given by 
Ea + AEa- The closest states generally have the great- 
est effect, and therefore most terms in the sum can be 
ignored. The T — 2 state in '^^S can be shifted by both 
T = and T = 1 states, while the states in ^^P and ^^Cl 
can only be shifted by T = 1 states (T = levels do 
not exist in these nuclei). Table II shows the significant 
contributing levels to energy shifts in the A ~ 32 quintet 
for all three interactions. 

The excitation energy of the T = 2 0+ state in 
32s is 11.885 (11.867,12.011) MeV with the isospin- 
nonconserving USDB (USDA,USD) interaction, in rea- 
sonable agreement with the experimental value of 12.048 
MeV. However, the nearest T = 0+ state is 182 (252) 
keV higher in energy for USDB (USD) but 33 keV lower 
in energy for USDA. There are no known experimental 
0+ states above the T = 2 state at 12.048 MeV; the 
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TABLE II: Energy shifts in the A = 32 quintet for the three 
interactions, where AEi,di, and are the contributions to 
the energy shift, d, and e coefficients, respectively, from state 
i (i.e., AEa ~ The algebraic solution to the five 

term IMME gives exact values of d and e for each interaction. 
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nearest known experimental 0^ level is a T = state 
118 keV below the T ~ 2 state. There is also an exper- 
imental O"*" state 180 keV below without an assigned T 
value which could correspond to the T — 1 state seen 
377 (302,195) keV below for USDB (USDA,USD). The 
proximity of these two states as calculated via the differ- 
ent empirical interactions has no deep underlying cause, 
but is rather an incidental effect due to the configura- 
tions of the states. Because the shift varies inversely with 
the energy, the proximity of the mixing state determines 
the size of the shift and the observed deviation from the 
three-term IMME. This energy difference (in conjunction 
with the size of the isospin-mixing matrix element) de- 
termines the size of the necessary d or e coefficient. In 
an algebraic solution of an isobaric quintet to the five 
term IMME (including both ^(r^)^ and e{T^)^ terms), 
the USDB and USD calculations, as well as the experi- 
mental data, result in a small ( < 0.1 keV) e coefficient. 
The USDA calculations require an e coefficient of 1.4 keV, 
even larger than the necessary d coefficient. Since the 
small energy difference in the USDA level scheme results 
in strong mixing in ^^S, in opposition to the experimental 
data, we will favor the USDB result. 

The decay from the T = 2 state in ^^Cl occurs via 
two processes: 7 decay and proton emission. The pri- 
mary channel of 7 decay is an Ml transition to the first 



excited l"*" state in ^^Cl with a branching ratio of 94%, 
using the USDB interaction and free gyromagnetic fac- 
tors. The calculated width is 1.11 eV, in comparison 
to the experimental value of — 1.8(5) eV. The isospin- 
forbidden proton transition decays to ^^S. The reaction 
has Q ~ 3.45 MeV for both USDB and experiment. Since 
the transition is from a J'^ = 0^ level, and the proton 
has j = 1/2+, 3/2+, 5/2+ in the sd shell, decay can only 
occur to levels in ^^S with those values of J. Five such 
levels have < Q, all with T = 1/2 using the isospin- 
conserving part of the Hamiltonians. With the inclusion 
of the CIB and CSB interactions, isospin mixing occurs 
in both the parent ■^^Cl and daughter '^-'^S nuclei. The 
small spectroscopic factors for the isospin-forbidden tran- 
sitions are shown in Table III. If isospin mixing is only 
included for ■^^Cl, the spectroscopic factors are larger. It 
is therefore important to include mixing in both nuclei 
to account for interference effects in the wavefunctions. 
The decay widths to states in '^^S were determined by 

Fp^^C^^F^p, (3) 

3 

where C^S are spectroscopic factors and Tgp is the single 
particle width of the resonance peak in the reaction '^^S 
-I- p — >■ '^^Cl. The single particle widths were calculated 
from scattering phase shifts in a Woods-Saxon potential 
|21) . with the potential depths chosen to reproduce 
the resonance energies. The results for the five levels are 
shown in Table III, but only the 1/2+ ground state and 
first excited state (3/2+) contribute to the decay. There- 
fore, Tp — 41.4 eV, in comparison to the experimental 
value of 20(5) eV from [S]. There is a large uncertainty 
in the 1/2+ level, the most important contribution, for 
two reasons: (i) the calculation of F^p was determined by 
doubling the width at half max on the low-energy side of 
the resonant peak, due to the large tail in the resonant 
reaction at high energy and (ii) the spectroscopic factor 
changes by a factor of four depending on the interaction 
used, with C^S values of 0.00003 (0.00008,0.00012) for 
the USDB (USDA,USD) interactions. Using the same 
single particle widths for the USDA interaction gives 
Fp = 92.2 eV. The isospin-mixing matrix element be- 
tween the T = 2 state in '^^Cl and the T = 1 state below 
it is approximately the same for all three interactions, 
as seen in Table II. The energy change difference varies 
from 262 to 439 keV, however, where the range is on the 
order of the 150 keV global rms energy deviation for the 
three interactions. However, a difference of 100 keV in 
the energy denominator can significantly affect the width 
of proton decay due to the amount of mixing. The USDB 
energy difference is greatest and has the smallest proton 
decay width, in better agreement with experiment. 

In the /?+ decay from the ground state of "^^Ar to 
the T = 2, J'^ = 0+ state in '^'^Gl (both members of 
the T = 2 multiplet), the ft value differs slightly from 
the expected value due to the isospin mixing in '^^Cl. 
The only significant contribution comes from the T=l, 
= 0+ state below, the same state infiuential in the 
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TABLE III: Widths of the isospin-forbidden proton decay for 
^^Cl using the USDB interaction. 
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isospin-forbidden proton decay to '^^S and in the devia- 
tion from the three-term IMME. The calculated value 
of 5c is 0.27% (0.40%,0.63%) for the USDB (USDA, 
USD) interactions, allowing us to quote a theoretical 
value of 0.43(20)% from an average of the three calcu- 
lations. Again, the energy difference results in a range of 
results for an isospin-mixing effect. From [6j, 5c should 
be a sum of a charge-dependent mixing contribution (cal- 
culated here) and a radial overlap component (1.4% from 
[5]). The sum of these two contributions, or 1.8%, agrees 
with the experimental values of 2.1(8)% [5] and, based 
on a new mass of ^^Cl, 1.8(8)% (TH]. In all three cal- 
culations, nearly the entire remaining strength feeds to 
the Oj state in ^^Cl, the T — 1 state shown in Table II. 
The transition to this state from the ground state of Ar 
might be accessible experimentally. 

The experimental data for the A = 32 multiplet differs 
from the IMME fit with three terms, requiring another 
term for an adequate fit. Using the new direct measure- 
ment of '^^Si [T?, a d coefficient of 0.93(12) is necessary. 
With the indirect mass of "^^Si [10], [T], the necessary d 
coefficient is 0.53(11). The three term IMME similarly 
does not reproduce the behavior of the masses using three 
different sd interactions. The calculated d coefficient is 
0.28 (0.30,0.39) for the USDB (USDA,USD) interactions. 
The USDA calculations also result in a large e coefficient 
due to the proximity of a T = state to the T — 2 state 



in ^^S. There is an inherent uncertainty in our method 
regarding the shift due to isospin mixing on account of 
the global rms deviation of 150 keV of empirical interac- 
tions. We gain information from using multiple interac- 
tions, but rely on experiment to constrain our choice of 
interaction for comparison. With the USDB interaction, 
the decay of the = — 1 state in the multiplet occurs 
primarily by proton emission with Tp = 41.4 eV, but the 
gamma decay with = 1.11 eV cannot be neglected. 
The proton decay width is approximately double the ex- 
perimental value, while the gamma decay width is in rela- 
tively good agreement with experiment. The theoretical 
Fp result varies significantly with the interaction used, 
suggesting a large uncertainty in the calculated value. 
Regardless of the interaction chosen, the mixing of the 
T = 1 and T = 2 states in •^^Cl causes a nonzero isospin- 
forbidden proton decay to T = 1/2 states in ^^S. The 
mixing of these same states also accounts for the devia- 
tion in the ft value of the /3+ decay of the ground state of 
•^^Ar. The isospin-breaking correction 5c = 1.8% agrees 
with the experimental value. 

These observed aspects of isospin mixing occur in re- 
lation to the proximity of the levels, separate from the 
correlation between the configurations of the states. In 
the event of small energy differences and non-negligible 
isospin-mixing matrix elements, the effects described 
above will be seen. The commonness of the fulfillment of 
these two requirements in other multiplets, such that ef- 
fects of isospin mixing occur, cannot be determined with- 
out more accurate theoretical energies or more complete 
experimental level schemes. 
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